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The physical basis for observed long-lived electronic coherence in photosynthetic light-harvesting 
systems is identified using an analytically soluble model. Three physical features are found to be 
responsible for their long coherence lifetimes: i) the small energy gap between excitonic states, 
ii) the small ratio of the energy gap to the coupling between excitonic states, and iii) the fact 
that the molecular characteristics place the system in an effective low temperature regime, even 
at ambient conditions. Using this approach, we obtain decoherence times for a dimer model with 
FMO parameters of ~ 160 fs at 77 K and ~ 80 fs at 277 K. As such, significant oscillations are 
found to persist for 600 fs and 300 fs, respectively, in accord with the experiment and with previous 
computations. Similar good agreement is found for PC645 at room temperature, with oscillations 
persisting for 400 fs. The analytic expressions obtained provide direct insight into the parameter 
dependence of the decoherence time scales. 



Electronic energy transfer is ubiquitous in nature and 
its dynamics and manipulation is of special interest in di- 
verse fields of physics, chemistry, biology and engineering. 
Under natural conditions loss of coherence is expected to 
occur on ultrashort times scale due to interaction with 
the environment. For example, results on betaine dye 
molecules [1] and on femtosecond dynamics and laser con- 
trol of charge transport in trans-polyacetylene [2] suggest 
that these time scales are very short, ^--^2.5 fs and ~3.7 fs, 
respectively. At high temperatures and for weak cou- 
pling to the environment, a classical treatment of the 
thermal fluctuations suggests that this time scale can be 
determined as tq = ^f9- j^Xk^T where A is the system 
reorganization energy [3, 4]. Based on this expression, 
the dephasing time for photosynthetic complexes, the sys- 
tems of interest in this paper, (wherein a typical value of 
the reorganization energy is A = 130 cm~^), can be esti- 
mated [5] to be TQ = 45 fs at T = 77 K and tq = 23 fs 
and T = 294 K. 

By contrast, recent experiments in photosynthetic com- 
plexes such as the FMO complex [6, 7] and the PC645 
complex [8] have found that electronic coherences among 
different chromophores survive up to 800 fs at 77 K [6] 
and up 400 fs at room temperature [7, 8]. This surprising 
observation and its possible consequences for biological 
processes have been discussed extensively [5-16] and very 
elaborate models have been developed in order to under- 
stand the underlying dynamics [11, 13-16]. Interestingly, 
despite the diversity of approaches and techniques, most 
[5, 8, 10, 11, 13-16] now predict long-lived coherences 
on the same times scales as found experimentally [6-8]. 
This suggests that the underlying physical features arc 
correctly contained in these approaches. However, the 
sheer complexity of these computations has limited one 
from identifying these essential physical features. 

In this paper we present a simple analytic approach 
that provides deep insights into the long lived coherences 
in the evolution of FMO complexes [5-8, 10, 11, 13, 15] 



and PC645 [8, 14] and allows for the identification of 
central characteristics responsible for these long lived co- 
herences. Our analysis identifies the small effective tem- 
perature of the system (see below), the very small, but 
nonzero, energy gap between exciton states, and their 
coupling as the basic elements behind these long lifetimes. 
Given these conditions, we show that the lifetimes are not 
"surprisingly long" . As such, the challenge now reverts to, 
for example, obtaining an atomistic understanding [15] of 
the origins of these parameter values. 

Consider first results for the Fcnna-Matthews-Olson 
(FMO) complexes, in particular the FMO pigment- 
protein complex from Chlorobium tepidum [6, 7]. This 
is a trimer consisting of identical, weakly interacting 
monomers [17]. Each weakly interacting FMO monomer 
contains seven coupled bacteriochlorophyll-a (BChla) 
chromophores arranged asymmetrically, yielding seven 
nondegenerate, delocalized molecular excited states (ex- 
citons) [6, 7]. Since the electronic coupling between the 
BChla 1 and BChla 2 is relatively strong in comparison 
with the other coupling strengths [17], we can approxi- 
mate the dynamics of the excitation as given by a dimer 
composed by BChla 1 and BChla 2. This, often adopted, 
approximation is utilized below. 

We consider the dimer (see 1) to be described by the 
Hamiltonian [18, 19] 
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is the reaction field op- 



where Ri = J2a ^a, 

erator for molecule i, Bi = ^^ huja,iO-a i^a,i is the energy 
stored in the solvent cage of molecule i and Sjij is the dif- 
ference between the dipolc moment of the chromophore 
i in the ground and excited states [18, 19]. The first 
two terms in 1 are the contributions from the Individ- 
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FIG. 1. Left hand side: The pair of interacting chromophores. 
Right hand side: the effective light harvesting two-level sys- 
tem formed from the pair of interacting chromophores [19, 20]. 



ual sites and the third term is the A coupling between 
them. The subsequent terms describe the system-bath 
coupling. Following references 19 and 20, the Hamil- 
tonian in 1 can be written with respect to the basis 

{|5i) ® I52), I51) ® 1^2), |ei) ® 1^2), |ei} ® |e2}} describing 
the state of the two chromophores, i.e. 
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where e± = ei ± £2, and V± = SfiiRi ± (5^2 ^2- 

Since under excitation by weak light only the singly 
excited states need to be taken into account, we can iden- 
tify [19, 20] the active environment coupled 2D-subspace 
as {|ei) ® 1^2), I51) ® |e2)}. In this central subspace of 
2, the effective interacting biomolecular two-level system 
Hamiltonian reads 



H = ( y CT^ + fiAa^ 
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y^ huja,ia'l^^aa,i (3) 
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where e = e_ and V = V^. This is schematically illus- 
trated in 1, where A is the associated "tunneling energy", 
between the new basis states |ei) (g) |g2) and \gi) ®\e2)- 

Given the biophysical nanostructure composition [19], 
we can assume that the two bath are uncorrelated 
[aQ,i,aI,2] = 0) which implies that {Ri{t")R2(t')) = 
{R2{t")Ri{t')) = 0. Hence, 3 can be written in the stan- 
dard form of the spin-boson model [19] 



H 



(4) 
where the hp includes harmonic oscillators coupled to 
both chromophores, with couplings gp. 

The environment is fully characterized by the spec- 
tral density J(w) = '^^d'a^i^ ~ ^ci)^ being a quasi- 
continuous function for typical condensed phase applica- 
tions that determines all bath-correlations that are rel- 
evant for the system [21, 22]. For Ohmic dissipation 



J{u)) ~ 2Kiuc~'^/^'^ ^ where the dimensionlcss parameter 
K describes the damping strength and uiq is the cut-off 
frequency [21, 22]. An Ohmic spectral density is a useful 
choice for, e.g. electron transfer dynamics or biomolecu- 
lar complexes [19, 23]. The parameter K is related to the 
reorganization energy A by means of A = 2KhuJc and the 
phonon relaxation time is given by r = tt/{2ujc) [11, 21]. 

The Hamiltonian in 4 has been extensively studied in 
literature (cf. Chaps. 18-22 in Ref. 21 and references 
therein). The parameter range within which the light 
harvesting systems of interest lie allows for the use of 
the non-Markovian non-interacting blip approximation 
(NIBA) plus first order corrections in the interblip cor- 
relation strength, i.e. an enhanced NIBA approximation. 
This approximation is valid for weak system-bath cou- 
pling and for e/2A < 1, over the whole range of tempera- 
tures (see Chap. 21 in Ref. 21), and provides simple and 
accurate analytic expressions for relaxation and decoher- 
ence rates. 

In the case of FMO, the energy gap e is (315 — 
240) cm~^ = 75 cm"-'^ while the coupling energy A corre- 
sponds to 87.7 cm~^ (see [11, 16] and references therein). 
In accord with references 11 and 16, the reorganization 
energy in this case is A w 35 cm~^ and the phonon 
relaxation time is t = 50 fs. Hence, for this case we 
get K ^ X/(2nuic) = XT/{hn) = 0.105. Additionally, 
at T = 77 K, 2A/kBT « 3.28 while at T = 277 K, 
2A/fcBT w 0.911. I summarizes the parameters of the 
present analysis, where we have used the same wc value 
at both temperatures. Clearly these parameters place the 
system within the domain of accuracy of the enhanced 
NIBA approach. We emphasize that this selection of pa- 
rameters is widely used in the literature [11, 16], and it is 
not chosen to fit our model to the experimental results. 
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TABLE I. Parameters used for dimer formed of BChla 1 and 
BChla 2 at T = 77 K (first row) and T = 277 K (second row). 



The high temperature limit in this approach is given 
by temperatures well in excess of Tb = fi.(A^ff-|-e^)^/^/fcB, 
where 

Aeff = [r(l - 2i^)cos(7ri^)]i/2(i-^')(AK)^/(i-^)A, 

where in our case A — 2A. 

For the set of parameters listed in I, we find Tb ~ 
288 K. Hence the FMO experiments, at 77 K and 277 
K, are in the low temperature regime, T < Tb. In this 
regime, the Rabi frequency fl, the relaxation rate 7r and 



the decoherence rate 7 are given by [21] 

r!2 = A2ff[l-2SRu(iAb)]+e' (5) 

^ A2 

f^(Ab) (6) 
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respectively, where Ab ~ V'^la + ^^ ^'^"^ ^(•^) ^ 
i/o°°dw;;j^(coth(fej/2fcBT) - 1) while S{uj) = 
J(a;) coth(?jw/2fcBT) is the noise power. Non-Markovian 
corrections are aheady included in 5 and 7. For the par- 
ticular case ol Ohmic dissipation adopted here [21], 

(8) 
+ 2KAls [»V(i/iAb/27rfcBT) - ln(nAb/27rfcBr)] 

7,. = 7:K coth(SAb/2A;BT) A^fj/Ab, (9) 

7 = 7r/2 + 2TTK{eyAl)kBT/h, (10) 

where tpi^) is the digamma function. 9-10 provide sim- 
ple analytic expressions for the desired rates. 

With this set of expressions, and the parameters given 
in 1, we find 27rf}-i = 163 fs, 7-1 = 90 fs, 7-1 = 153 fs at 
T = 77 K, while 27rf7-i = 151 fs, 7^1 = 45 fs, 7-^ = 69 fs 
at T = 277 K. Despite the simplicity of the model, the 
resultant dynamics [given analytically, for general spec- 
tral densities, in Eqs. (21.79), (21.170) and (21.173) of 
Ref. 21] is depicted in 2 and describes the survival of 
coherences on the correct time scale and in good agree- 
ment with recent results[ll, 15, 16]. In 2, only the decay 
of the excitations is absent, since coupling to other chro- 
mophores [11, 15, 16] is here neglected. The global decay 
in 2 is a result of two processes: one associated with the 
71- relaxation of the dynamics, and a second one related 
to the loss of coherence due to the presence of the thermal 
bath, associated with 7. 

To examine the parameter dependence, 3 shows the re- 
laxation rate 71- and the decoherence rate 7 as a function 
of the energy splitting, e, for fixed A, using the param- 
eters displayed in I. Clearly in this "low temperature" 
regime, i.e., T < T/,, the larger the ratio e/2A the shorter 
the dephasing time and concomitantly the longer the re- 
laxation time. These time scales are also compared in 3 
with the ones predicted by 7^ ~ 2'K{kBT / h)\ / hwc used 
in references 4, 7, 10 and with 7g = I/tq discussed ear- 
lier, both of which are clearly far too short. 

As a second example we consider results on marine 
algae [5], in particular the results for PC645, which has 
been studied experimentally [8] , and numerically in great 
detail [14]. PC645 contains eight bilin molecules cova- 
lently bound to the protein scaffold. A dihydrobiliverdin 
(DBV) dimer is located at the center of the complex 
and two mcsobiliverdin (MBV) molecules located near 
the protein periphery give rise to the upper half of the 
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FIG. 2. Time evolution of the population of the electronic 
manifolds jei)(g)|(;2)(ff2|'8(ei| (blue curves) and |(;i)(g)|e2)(e2|® 
{gi I (red curves) for the FMO dimer formed of BChla 1 and 
BChla 2 using the parameters displayed in I. Continuous 
curves depict the evolution at T = 77 K whereas dashed 
curves are at T = 277 K. We assume that at t = 0, the 
full excitation is localized in |ei) ® |52)(<?2| ® (ei|. 
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FIG. 3. Relaxation time 7^ ^ (red line) and decoherence time 
7~^ (continuous blue line) as functions of e/2A, for fixed 
A = 35 cm"^ and at T = 77 K. The dashed blue line, 
70 — 2Tr {ksT /h)X/ hue , denotes the dephasing rate used in 
references 4, 7, 10 and th e dotdashe d blue line denotes the 
decoherence time ra = \/f^ l2\k^T used in references 3-5. 
Fixed parameters as in I. 



complex absorption spectrum. Excitation of this dimer 
initiates the light harvesting process [8]. The electronic 
coupling between the closely spaced DBVc and DBVd 
molecules is ^^320 cm~^ and this relatively strong cou- 
pling results in delocalization of the excitation, yielding 
the dimer electronic excited states labelled DBV-I- and 
DBV—. Excitation energy absorbed by the dimer fiows 
to the MBV molecules which are each 23 A from the clos- 
est DBV, and ultimately to four phycocyanobilins (PCB) 
that absorb in the lower-energy half of the absorption 
spectrum [8]. 




FIG. 4. Time evolution of the population of the electronic 
manifolds |ei)8|(;2)(52|® (ei| (blue curve) and |gi)®|e2)(e2|® 
{gi \ (red curve) for the PC645 DVB-dimer. 



The exciton states related to DBVc and DBVd arc 
mainly composed of DBV— and DBV+ that arc antisym- 
metric and symmetric linear combinations of the DBV 
sites, though they also contain small contributions from 
the other bilin sites [14]. This allows us to concentrate 
only on a dimer, as in the previous example, here formed 
of DBVc and DBVd. In this case, a Debye-Ohmic spec- 
tral density, J(w) = 2Xujt/{1+uj'^t^), is more appropriate 
[14]. 

For the chromophores DVBc and DVBd the energy 
gap is (17116 - 17034) cm"! = 82 cm^^ with the cou- 
pling energy corresponding to 319.4 cm~^ [8, 14]. In 
accord with Ref. 14, the reorganization energy A in this 
case is ~130 cm^^, with the shorter of two relaxation 
times being r = 50 fs. Additionally, at T == 294 K, 
the temperature at which the experiment was performed, 
2A/kBT = 3.1262. For this set of parameters, we have 
Tb ~ 926 K, a high temperature that is consistent with 
the high frequencies involved in the present case [24]. 
Hence the PC645 experiment, at room temperature, is 
also in the effective low temperature T < Tb domain. 
For the case of the Debyc-Ohmic spectral density we 
have (via a relatively simple numerical computation) that 
27rf}-i =: 49 fs, 7-1 =: 76 fs, 7"^ = 88 fs with the asso- 
ciated evolution shown in 4. These results are in very 
good qualitative agreement with both the long lived ex- 
perimental coherence time scales and with the recent in- 
tricate solution to the master equations [14]. 

Considering that the long time scales emerge natu- 
rally here from the system parameters, why were far 
shorter decoherence time scales originally expected for 
these systems? To see this, note that in molecular sys- 
tems, dynamics is often studied between different elec- 
tronic eigenstates of the system, separated by greater 
than ^ 10^ cni^^, with no coupling between them. In 
such cases, the dephasing time from 9 and 10 would be 
extremely short. By contrast, in the case of photosyn- 



thetic complexes, energy transfer occurs between exci- 
ton states that are close in energy and additionally are 
coupled. This generates a small value for the ratio e/A 
which in turn is responsible for longer dephasing times 
(see 3). Additionally, expression such as tq, which are 
often used to estimate rates, are only valid at high tem- 
perature, k^T ^ huJc, and at short times, t < ojq . Un- 
der conditions when the expression for tq is valid, the 
bath modes can be treated classically [25], as in refer- 
ences 1 and 2. When one is not in the appropriate regime, 
the classical evolution of the bath underestimates quan- 
tum coherence effects [25] because at low temperatures 
quantum fluctuations overcome thermal fluctuations [3]. 
Hence, estimates based on tq are unreliable. Similarly, 
70 ~ 2tt {k^T / h)\ / tiwc estimates also provide an inade- 
quate representation of the true physics and associated 
dependence on system and bath parameters, and result 
in decoherence times that are severely underestimated 
and misleading. 

In summary, a proper spin-boson treatment of elec- 
tronic energy transfer in model photosynthetic light har- 
vesting systems has been shown to give analytic results 
with long coherence lifetimes that are in very good agree- 
ment with experiment and with other, far far more com- 
plex, computations. The analytic form allows an analy- 
sis of the parameter dependence of the decoherence times, 
and shows that the observed long lifetimes arise naturally 
in the effective low temperature regime and for appropri- 
ate ratios of the energy splitting to the coupling strength 
and are, in this sense, not surprisingly long. Further, 
the model has predictive power, allowing one to identify 
other parameter ranges over which longlived coherences 
will exist. 

A more detailed analysis will involved extending this 
model to the 7 or 8-site system and to second order cor- 
rection in the intcrblip approximation; work in this di- 
rection is in progress. Finally, in addition to examining 
other light harvesting systems, the possible application of 
this approach to other systems displaying long-lived co- 
herences, such as superpositions of two excitons in quan- 
tum dots [26] is of interest. 
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